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The domain model for the QCD vacuum has previously been developed and shown to exhibit 
confinement of quarks and strong correlation of the local chirality of quark modes and duality of 
the background domain-like gluon field. Quark fluctuations satisfy a chirality violating boundary 
' conditions parametrized by a random chiral angle ctj on the j — th domain. The free energy of 

| an ensemble of N — > oo domains depends on {ctj , j = 1 . . . N} through the logarithm of the quark 

determinant. Its parity odd part is given by the axial anomaly. The anomaly contribution to the free 
£\j ■ energy suppresses continuous axial (7(1) degeneracy in the ground state, leaving only a residual axial 

Z(2) symmetry. This discrete symmetry and flavour SU(Nf)L X SU(Nf)n chiral symmetry in turn 
are spontaneously broken with a quark condensate arising due to the asymmetry of the spectrum of 
Dirac operator. In order to illustrate the splitting between the r/ from octet pseudoscalar mesons 
realised in the domain model, we estimate the masses of light pseudoscalar and vector mesons. 
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I. INTRODUCTION 

m ■ „ 

A mechanism which simultaneously provides for confinement of colour, spontaneously broken chiral symmetry and 
t-H | a resolution of the Ua(1) problem remains one of the open problems in QCD today. Partial solutions [1] based on 
■ specific semi-classical or topologically stable configurations can go someway to manifesting this triplet of phenomena 
but either founder on generating all three or in allowing for an effective model of the vacuum from which hadron 
spectroscopy can be derived. In any case, one expects that topological objects of various dimensions - point-like, 
string-like, and sheet-like - should contribute [2] complete with significant quantum fluctuations in a way that would 
Oh, be difficult to describe via an interacting microscopic model. In this paper, we continue the exploration of the "domain 
model" for the vacuum, originally proposed in [3] , as a scenario for simultaneous appearance of all three phenomena: 
| confinement, spontaneous chiral symmetry breaking via the appearance of a quark condensate and a continuous 
' SU(Nf)i, x SU(Nf)n degeneracy of the vacuum for Nf massless quarks, but without a J7a(1) continuous degeneracy 
| of ground states that would be indicative of an unwanted Goldstone boson. The purpose of studying a model of this 
type is to identify the typical features of relevant nonperturbative gluonic configurations. Such configurations would 
provide for as many gross features of nonperturbative QCD as possible. But this should preserve simultaneously the 
well-studied short distance regime and the model should be expressed in terms of quark-gluon degrees of freedom as 
well as in terms of colourless hadron bound states. 

The model under consideration provides for confinement of both static (area law) and dynamical (propagators are 
entire functions of momentum) quarks [3]. It also displays specific chiral properties of quark eigenmodes. Namely, 
as will be discussed in more detail below, the spectrum of the Dirac operator is asymmetric with respect to A — * — A 
and zero quark modes are absent, but the local chirality of all nonzero modes at the centre of domains is correlated 
with the duality of the background field [4] . This has been observed on the lattice [5] and is usually considered as an 
indication of spontaneous breakdown of flavour chiral symmetry. The purpose of this article is to study the details of 
chiral symmetry realization in the domain model. The nonzero quark condensate and axial anomaly are generated as 
a result of spectral asymmetry and definite mean chirality of eigenmodes. We compute the quark condensate, study 
the degeneracies of the minima of the free energy of the domain ensemble with respect to chiral transformations and 
estimate the spectrum of pseudoscalar mesons. 

The model is defined by a partition function describing an ensemble of hypersphcrical domains, each characterised 
by a background covariantly constant self-dual or anti-self-dual gluon field of random orientation. Summing over all 



* akalloni@physics.adclaidc.edu.au 
t ncdelko@thsunl.jinr.ru 



2 



orientations and both self-dual and anti-self-dual fields guarantees Lorentz and CP invariance. Quarks are confined as 
demonstrated in the original work [3]. On the boundaries of each hypersphere, fermion fluctuations satisfy a chirality 
violating boundary condition 

ifl(x)e ia ^il){x) = 4>(x) (f) 

which is 2n periodic in the chiral angle a. Here 77^ is a unit radial vector at the boundary. Integrating over all such 
chiral angles guarantees chiral invariance of the ensemble. As a consequence of Eq. (1) the spectrum of eigenvalues 
A of the Dirac operator in a single domain is asymmetric under A — > —A. Such asymmetries have been studied in 
other contexts, for example by [6]. In the case of the domain model, the above boundary conditions are combined 
with the (anti- )self-dual gluon field which leads to a strong correlation between the local chirality of quark modes at 
the centres of domains with the duality of the background gluon field [4]. In this paper we study how these aspects 
contribute to quark condensate formation and the pattern of chiral symmetry breaking. 

The vacua of the quantum problem associated with an ensemble of domains are the minima of the free energy 
determined from the partition function. The problem of the quark contributions to the free energy requires calculation 
of the determinant of the Dirac operator in the presence of chirality violating boundary conditions. For a choice of 
boundary condition with a — > —id — 7r/2 this problem has been tackled in [7] without a taking into account of the 
spectral asymmetry where the parity odd part of the logarithm of the determinant was identified as lndet(i P) ~ 2q$ 
with q the topological charge (not necessarily integer) of the underlying gluon field, namely the axial anomaly. 

Our first goal in this paper is to address the analogous problem for the specific gluon field relevant to the domain 
model taking into account the asymmetry of the spectrum. For the parity odd part we obtain 

lndet(i,0) ~ 2iq(a mod tt). (2) 

This result is consistent with [7] up to a contribution coming from the asymmetry spectral function. However, we 
obtain an additional parity even part which also turns out to be a-dependent. We consider this to be more an artifact 
of the incompleteness of our calculation than an established property of the determinant. 

In the partition function all possible sets of chiral angles {ot\, . . . , ajv} are summed ensuring the chiral invariance 
of the ensemble. Summation over all degrees of freedom besides chiral angles defines the free energy as a function of 
these chiral angles. In the limit N — > 00 the minima of the free energy density in {a\, . . . , a^v} determine the preferred 
chiral angles. More specifically, when self-dual and anti-self-dual configurations are summed, the anomaly Eq. (2) 
leads to a contribution to the free energy of the form — lncos(2<7arctan(tanai)) which vanishes when a — 0, n. The 
minima are degenerate with respect to discrete Z 2 chiral transformations. Each of these minima are characterised by 
a quark condensate of opposite sign, which arises due to the the spectral asymmetry. An infinitcsimally small quark 
mass removes the degeneracy between the two discrete minima, and a nonzero quark condensate is generated with 
value 

(ip(x)tp(x)} = -(237.8 McV) 3 

with no additional modifications of the two model parameters after fixing in the gluonic sector of the theory. This 
gives a model with the chiral Z2 discrete subgroup of 17a (1) being spontaneously broken, and not the continuous 
[/^(l) itself. In the absence of the mass term the ensemble average of i/jip correctly vanishes. A similar argument 
based on minimisation of the free energy and thereby a relaxation of the effective ^-parameter of QCD to zero is 
discussed in detail in [8] in the context of the strong CP-problem. 

Moreover, the form of Eq.(2) means that the free energy does not depend on flavour nonsinglet chiral angles when 
more than one massless quark flavours are introduced. This allows for the correct degeneracy of vacua with respect to 
continuous SU(Nf)i x SU(Nf)n chiral transformations. This vacuum structure implies the existence of Goldstone 
bosons in the flavour nonsinglet pseudoscalar channel but not in the singlet channel. To unveil more explicitly the 
singlet-octet splitting we analyse the structure of pseudoscalar correlation functions in the context of the domain 
model and estimate the masses of light pseudoscalar and vector mesons. The qualitative conclusion of this analysis 
is that the area law (confinement of static quarks) and the singlet-octet splitting in the model have the same origin: 
the finite range correlations of the background gluon field. 

In the next section we briefly review the model, followed by a summary of the properties of the spectrum of the 
Dirac operator in the domain-like gluon field. We then discuss in detail the calculation of the logarithm of the quark 
determinant for one massless quark flavour, including the role of spectral asymmetry in domains in giving the anomaly 
for the parity odd part. This is followed by an analysis of the symmetries of the ground state of the domain ensemble 
and the computation of the condensate. In section V we generalise the result to Nf massless flavours in order to 
verify the spontaneous breaking of SU(Nf)L x SU(Nf)n in the ensemble. The last section is devoted to calculation 
of meson masses. Details of calculations are relegated to the appendix. 
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II. THE DOMAIN MODEL 



For motivation and a detailed description of the model we refer the reader to [3] . The essential definition of the 
model is given in terms of the following partition function for N — > oo domains of radius R 



N 



Z=M lim WfdaJ 2ty«XtyM f VQ 1 S[D(B^)Q^]Afp[B {i \ Q(0] e -^° D [O (i) +B«^«,*«] (3) 

V,^ooA±7 ^ J 

where the functional spaces of integration Tq and T" 1 ^ are specified by the boundary conditions (x — Zi) 2 = R 2 

tUQ®(x)=0, (4) 
iff l (x)e la '^^ l) (x) = ip^(x), (5) 
^e ia ^Hffi(x) = -$®(x). (6) 

Here rn — nft a with the generators t a of SU C (3) in the adjoint representation and the Qj are chiral angles associated 
with the boundary condition Eq.(5) with different values randomly assigned to domains. We shall discuss this 
constraint in detail in later sections. The thermodynamic limit assumes V,N —> oo but with the density w _1 = N/V 
taken fixed and finite. The partition function is formulated in a background field gauge with respect to the domain 
mean field, which is approximated inside and on the boundaries of the domains by a covariantly constant (anti-)self- 
dual gluon field with the field-strength tensor of the form 

N 

n(i) hO) — r2 c 

Here Zj are the positions of the centres of domains in Euclidean space. 
The measure of integration over parameters characterising domains is 



2tt 



3.4,5 



/ d<Jl " ' = 4^r2 / ~V~ I dai I dtft f Mt Sin ° l f d ^ X/ S (£i~ — "t 1 ' >7r ) / ^ S ( Ui ~ 7rfc ) • • • ' ^ 

J n JV JO J J0 Jo ;=o,l,2 ' ^° fe=0,l 

where (6i,(fi) are the spherical angles of the chromomagnetic field, Wj is the angle between chromoelectric and 
chromomagnetic fields and & is an angle parametrising the colour orientation. 

This partition function describes a statistical system of the domain-like structures, of density v^ 1 where the volume 
of a domain is v = ir 2 R i /2. Each domain is characterised by a set of internal parameters and whose internal 
dynamics are represented by fluctuation fields. Most of the symmetries of the QCD Lagrangian are respected, since 
the statistical ensemble is invariant under space-time and colour gauge transformations. For the same reason, if the 
quarks are massless then the chiral invariance is respected. The model involves only two free parameters: the mean 
field strength B and the mean domain radius R. These dimensionful parameters break the scale invariance present 
originally in the QCD Lagrangian. In principle, they should be related to the trace anomaly of the energy-momentum 
tensor [9, 10] and, eventually, to the fundamental scale Aqcd- Knowledge of the full quantum effective action of QCD 
would be required for establishing a relation of this kind. 

A straightforward application of Eq. (3) to the vacuum expectation value of a product of n field strength tensors, 
each of the form 



JV 



j 

gives for the connected n-point correlation function 

N 

(B;i Vl ( Xl )...B; :Vn {x n )) = J J?L Jd* j nV^...n^B$ Vl ...B$ Vn 



o 

x 9(1 - ( Xl - zjf/R 2 ) ...9(1- (x n - z f/R 2 ) 

— gny.oi...o„ w I \ 
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where the tensor t is given by the integral 

f ai...a n — / A Q>i n (j)a n u(j) u(j) 

and can be calculated explicitly using the measure, Eq. (7). This tensor vanishes for odd n. In particular, the integral 
over spatial directions is defined by the generating formula 



sin J2B 2 [J fll/ J tll/ ± J^J^} 



— / dpj / '/".MU'V'''.. — ■ 



The translation-invariant function 



~ n (x u x n ) = 1 J d 4 z6{l - ( Xl - z) 2 /R 2 ) ... 9(1 - (x n - zf/R 2 ) (8) 

can be seen as the volume of the region of overlap of n hyperspheres of radius R and centres (x\, . . . ,x n ), normalised 
to the volume of a single hypcrsphere v = ir 2 R 4 /2, 

S„ = 1, for xi = ■ ■ ■ = x n . 

It is obvious from this geometrical interpretation that S„ is a continuous function and vanishes if the distance between 
any two points \xi — Xj\ > 2R; correlations in the background field have finite range 2R. The Fourier transform of 5 n 
is then an entire analytical function and thus correlations do not have a particle interpretation. It should be stressed 
that the statistical ensemble of background fields is not Gaussian since all connected correlators are independent from 
each other and cannot be reduced to the two-point correlations. 

Within this framework the gluon condensate to lowest order in fluctuations is AB 2 , the absolute value of the 
topological charge per domain reads q = B 2 R 4 / 16 and the topological susceptibility turns out to be x = B 4 R 4 / 128n 2 . 
An area law is obtained for static quarks. Computation of the Wilson loop for a circular contour of a large radius 
L >• R gives a string tension a = Bj{nBR 2 ) where / is given for colour SU{2) and SU(3) in [3]. The area law 
emerges due to the finite range of background field correlators Eq. (48). On the other hand, the model cannot account 
for such a subtle feature as Casimir scaling: the adjoint Wilson loop naturally shows perimeter law, but trivially 
because of the abelian character of the domain mean field. 

Estimations of the values of these quantities are known from lattice calculation or phenomenological approaches and 
can be used to fit B and R. As described in [3] these parameters are fixed to be VB = 947MeV, R = (760MCV)- 1 = 
0.26fm with the average absolute value of topological charge per domain turning out to be q ~ 0.15 and the density of 
domains v -1 = 42fm~ 4 . The topological susceptibility is then \ = (197MeV) 4 , comparable to the Witten-Veneziano 
value [11]. This fixing of the parameters of the model remains unchanged in this investigation of the quark sector. 
The quark condensate at the origin of a domain where angular dependence drops out was estimated in paper [3] with 
a result of -(228MeV) 3 . 



The eigenvalue problem 



III. DIRAC OPERATOR AND SPECTRUM 



Pil>{x) = XiP(x), 
ifl{x)e lai5 4>{x) = i/>(x), x 2 = R 2 



was studied in [4]. Dirac matrices are in an anti-hermitean representation. For a assumed to be real a bi-orthogonal 
basis has to be constructed. Solutions can be labelled via the Casimirs and eigenvalues 

K 2 = K^^ + l), fc = 0,l,...,oo 

K l,2 -> m l,2, 

mi, 2 = -fc/2,-fc/2+l,...,ft/2-l,fc/2, 
corresponding to the angular momentum operators 

Ki i2 = l(L±M) 
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with L the usual three-dimensional angular momentum operator and M the Euclidean version of the boost operator. 
The solutions for the self-dual background field are then 

feu = * rfXkrnr + VkZr > ( 9 ) 

where \ an d <P must both have negative chirality in the self-dual field and k is related to the polarisation of the field 
defined via the projector 

K = 7V+£ K + 7V_£_ K (10) 

with 

N ± = \{l±h/\h\\ E± = i(l±SB/B) 

being respectively separate projectors for colour and spin polarizations. Significantly, the negative chirality for \ 
and ip is the only choice for which the boundary condition Eq.(5) can be implemented for the self-dual background. 
The explicit form of the spinors \ an d <p can be found in [4], where it is demonstrated that the eigenspinor Eq. (9) 
has definite chirality at the centre of domain correlated with the duality of the gluon field. The boundary condition 
reduces to 

X = -e T "V, X = Ve™, ^ = R\ (11) 

where upper (lower) signs correspond to ip and x with chirality which, using the solutions, amounts to equations 
for the two possible polarisations, for A^ + : 



e~ ia M (k + 2 — A 2 , k + 2, z ) - ^ 
v ' lA 



and for A fc " 



2 



M (k + 2 - A 2 , k + 2, z ) - fc + 2 A M (k + 3 - A 2 , k + 3, z ) 



= 0, (12) 



e~ %a M (-A 2 , k + 2, zo) + ^x?M (l - A 2 , k + 3, z Q ) = (13) 

where zq = BR 2 /2 and A = A/ \J 2B. For the present work Eqs.(12,13) are the starting point, from which we see 
by inspection that a discrete spectrum of complex eigenvalues emerges for which there is no symmetry of the form 
A — > — A. For given chirality and polarisation and angular momentum k, an infinite set of discrete A are obtained 
labelled by a "principal quantum number" n. 

IV. QUARK DETERMINANT AND FREE ENERGY FOR A SINGLE DOMAIN 

A. Massless case 

We consider the one-loop contribution of the quarks to the free energy density F(B, R\a) of a single (anti-)self-dual 
domain of volume v = tt 2 R 4 /2 

exp{-vF(B,R\a)} = det a ( 



n 



AL(0) 



K,fc,n,mi 

= cxp{-C'(.s)} s=0 . (14) 

The normalization is chosen such that lim^^o F{B, R\a) = 0. The free energy is then F — w _1 C'(0). 

In the zeta-regularised determinant an arbitrary scale \i appears and it is convenient to work with scaled variables 



(3 = v2B/fi, p = ^R, i = \{B)/n, ^ = A(0)//i, 

and where the dimensionless quantity z = BR 2 /2 — /3 2 p 2 /A appears prominently Moreover it is convenient to 
analytically continue a — ► —I'd — n/2 to guarantee a real spectrum of eigenvalues of the Dirac operator for all We 
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shall regard the background domain field as being self-dual in the following. Then, Eqs.(12,13) can be recast into a 
form determining the rescalcd eigenvalues £, namely 



M(k + 2-e/P 2 ,k + 2,z) + e^ 



k 4- 9 - F 2 I R 2 

M(k + 2 - e/P 2 , k + 2,z)- \ Y M ( fc + 3 - ^//3 2 , k + 3, z) 



(15) 



for £ + and 

M(-e/P 2 , k + l,z) + e *^i^M(l - e/P 2 , k + 2, z) = (16) 

for £ . The zeta function £(s) breaks up into two parts [6], respectively symmetric (S) and antisymmetric (AS) with 
respect to £ — > — £ 

C(*) = Cs(s) + Cas(s), 

with 

Cs(«) = ^(l + e^ s )^( S /2), (17) 

CasW = \(l-e^) v { s ). (18) 

The assumption behind these representations is that the spectrum of eigenvalues can be well ordered according to 
the magnitude of the eigenvalues, |£|, which is certainly the case with the solutions to Eqs.(15,16). However there 
is an ambiguity which can be fixed by specifying whether the smallest-in-magnitude eigenvalue is either positive or 
negative. This determines respectively the sign choice in Eqs.(17,18). 

Thus the key quantities devolve into the zeta function for the squared Dirac operator and the asymmetry function 
respectively: 

Ws)= ^ (k+1 Amw-~w) (19) 



k,n 

's§"(«„(B)) sgn(S n (0)) 



kn (B)\ s I4 K J0)I 5 



(20) 



It should be stressed that in the presence of bag-like boundary conditions (g & n d Cas do not have the meaning of parity 
conserving and parity violating terms since a parity transformation in terms of eigenvalues is given by — ► — £(— 1?) 
and both spectral functions contain parity conserving and violating terms. Thus the determinant for a given parameter 
$ is defined by 

C'(0)=(^ 2 (0)±^(0)T^r?(0)), (21) 

with the normalization chosen in Eq.(14) such that C'(0) vanishes as B — > 0. 

Spectral sums over quantum labels N can be computed using a representation of the sum as a contour integral of 
the logarithmic derivative of the function whose zeroes determine the spectrum, 



E- 

N ^ 



1 '^in/fc), (22) 



? N 2m J T e di 

see [12], where the zeroes of /(£) = are £ = ^ N and the contour is chosen such that all zeroes are enclosed. With real 
parameter the poles lie on the real axis, and there is no pole at the origin for any ■&. By deforming the contour and 
accounting for the vanishing of contributions to the integral at infinity, the expressions arising from Eq. (22) can be 
transformed into real integrals. The following representations for the two spectral functions are eventually obtained 



A 2s sin(7rs) ^ r°° fit d 

h Jo ^ : 



, (s) .,.« Et .-. r^ (MM) , (23) 
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with \& and $ being the sum of contributions from the two polarisations, taking the form 



*(M|«.0) = ^ In 



A*(M|*) 



-,2tS> 



i?2(fc,i|z) 



re=± 



A 2 (k,t) 
A K (fc,i|z) 



e 2 »B 2 (k,t) 
■ie*B K {k,t\z) 



A K {k,t\z)-ie<>B K (k,t\z) 



and where 



A_(fc,t|z) 
S_(fc,t|z) 
A+(M|«) 
B+(k,t\z) 

A(k,t) 
B(k,t) 



M 



k 2 t 2 p 2 



Az 



ktp 



k + l,z 
k 2 t 2 P 2 



M 1 + 



2(fc+l) 

fc 2 iV 



A/ 



4z 



4z 

, fc + 1, — z 



,k + 2,z 



2z 
ktp 

2 k k\ 



k 2 t 2 n 2 



4z 



-M 



fc 2 t 2 p 2 
47" 



(ktp) k 
2 k kl 
(Hp) 



h{ktp), 
k h+i(ktp). 



(24) 



The Bessel functions Ik emerge from the limit B — > with the normalization of the determinant as specified above. 
The next step is to expand the confluent hypergeometric functions in 1/k, similar to the Debye expansion of Bessel 
functions, for example 



M 



(UUL,k + i,z) = C(t P ,k)J2- 



M n (tp,z) 
k n 



(25) 



The form of the prefactors and the M n (x, z) functions are given in Appendix A for the various Kummer functions 
appearing in the A K and B K . 

Before proceeding with more detail, let us give an overview of the subsequent steps. The expansions Eq. (25) are 
first inserted in £(s) and rj(s). Then the integrals over t can be evaluated term by term in the series in 1/k. The order 
of summation over k and n are then exhanged. This is the most subtle step, which we discuss below. But it means 
that now the sums over k can be read off in terms of the Riemann zeta function. The resulting C,pi (s) then has the 
structure 



C (s) = af ^E^l V Cr(2s + n - l)f(z 2 \n) + 5( p ,(s), 

ITT i 



(26) 



n=0 



where the term 8C, denotes those potentially present contributions coming from interchange of the order of summations 
over n and k. A similar structure appears for the asymmetry spectral function as well. Then we are interested in 
the decomposition of the resulting expressions around s = 0. In this limit the first term in Eq.(26) provides only a 
contribution from the n = 2 term and can be calculated with relative ease since only the lowest coefficients Mi and 
Mi in Eq. (25) contribute. However, as occurs in even simple problems [13], the second term in Eq.(26) is much more 
difficult to compute. To achieve this one needs to know coefficients f(z 2 \n) as a function of continuous variable n. 
Below we will elucidate on the contribution from the first term alone, bearing in mind the necessity of a complete 
analysis. 

Now in more detail, using the expansions of the type Eq.(25) as given in Appendix A, we arrive at the following 
expansions of ^ and $ in 1/k: 



k 2 



+ 0(l/k 3 ), 
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with the coefficient functions 



2 y(y-l)(3y + l) 

3 (y + i) 2 ' 
y 2 (y-i) 



* 2 (y|0) 



[y + l-(y 2 -l) 2 e 2 "] 2 
x [(y 2 - l)(2y 3 + 6y 2 + 7y + 1) - 2(y + 1)(2 ? / + 5y 3 + by 2 -y+ l)e 2 " 
+(y - l)(2y 4 + 6y 3 + 5y 2 -Ay- l)e ii} ] , 
{l-yfe™ -{l + yf 



2ieVvT^y 2 - 



[l + y + (l-y) e 2 *] 



2 ' 



with y = 



, 1 . . We thus get for s << 1 

1+t 



C^ 2 («) 
X 2 



z 2 (l + 2slnp)s 



I 



00 dt_d_ 
t^di 

° dt_d_ 
W s Jt 



^i{y) = -s + 0(s 2 ) 
* 2 (y|0) = 2s 



ln2 + ln(l + e 2 ' 5 ) 



+ 0(s 2 ) 



and 



z 2 1 
V(s) = - — (l + s\np)-J 2 

ITT S 



J 2 = 2ie h 



i^2(y\$) = -4 k + 2Arctan(sinh(tf)] + 0(s 2 ) 
t s at 2 



upon performing the t— (or y— ) integrations and where we have used that 

Cr(0) = -^, Cr(1 + 2 S ) = 1+ 7 + 0( s ). 

The final results for the first term in Eq. (26) and its analogue in the asymmetry function rj(s) are then summarised 
in the following equations: 



^-ln2 + ln(l + e 2l >) 



z 2 z 2 

r)(0) = — H Arctan(sinh(tf)) 

2 7T 



and hence 



C'(0) = y 



4 2 



(27) 



+ In (1 + e 2 ®) =F iArctan(sinh(tf))] + ^'(0), 



Now we can straightforwardly continue $ — > zcn + 27r/2. The final result for the free energy F is complex with the 
imaginary part of the form 



Q?F = ±2gArctan(tan(a)) 



(28) 



where q is the absolute value of topological charge in a domain. This charge is not integer here in general but the 
anomalous term is wn periodic in a. This is the Abelian anomaly as observed within the context of bag-like boundary 
conditions by [7]. Its appearance here is in the spirit of the derivation by Fujikawa [14], where the phase appears as 
an extra contribution under a chiral transformation on the fermionic measure of integration. 
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However the analytic continuation of Eq.(27) also exposes an a dependent real part, 



F(a) = 



ln(l ± cos(a)) 



where the signs are according to the prescription in Eqs.(17,18). There is no reason to expect that anything other 
than the anomaly should appear in the logarithm of the determinant. Thus the likelihood is that this additional real 
part contribution should be cancelled by the contributions <5£ in Eq.(26), which remains to be verified. 

Below we assume that the anomaly Eq. (28) provides the entire result for the a-dependent part of the free energy 
of massless fermions in a domain, which leads to intriguing consequences. 

B. Massive fermions 

Inclusion of an infinitesimally small fermion mass leads to a modification of the free energy by a term which is 
linear in mass to leading order, namely [6] 



. m 



F = F m=0 + i—r](l). 



(29) 



Using again the representation Eq.(23) but now in the vicinity of s = 1 one can find the following representation for 
the summed contributions of both polarisations in the self-dual domain to the asymmetry function 



fc=i 



l + Af(l,fc + 2,z) 



k + 2 



M(l,k + 3,-z) 



(30) 



We next evaluate the asymptotic behaviour in k. A singular term as s — > 1 can be extracted, which turns out to be 
independent of the field (that is, B) and is canceled by the normalization. The final expression for the free energy 
density for a self-dual domain, including the contribution linear in mass, is 

F (sd) = 2i g Arctan(tan(a)) - e ia m -< ^ y, 

where we have used the suggestive notation 

1 A k 



E 



7T 2 i? 3 ^ k + 1 

fc=l,Z = Zi ,Zl ,Z2 



M(l,k + 2,z)- 



fc + 2 



M(l,k + 3,z) - 1 



(31) 



coming from 77(1) with the sum over z correponding to a colour trace. 

Let us summarise the results of this section as follows. The a dependent part of the free energy of a self-dual 
domain for massive quarks is complex with the following real and imaginary parts: 



F = $tF + i$sF 
$tF = —mcosa^,ipipy 

"^sF = 2-Arctan(tan(a)) — m sin a -< ipip >~ . 
The free energy of an anti-self-dual domain is obtained via complex conjugation. 

V. ENSEMBLE FREE ENERGY AND CHIRAL SYMMETRIES. 
A. One flavour case: quark condensate and C/a(1) 



(32) 
(33) 



Under the assumption that only the anomalous term depends on the chiral angle, the part of the free energy density 
T relevant for the present consideration of an ensemble of N — > 00 domains with both self-dual and anti-self-dual 
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configurations takes the form 



= A/"cxp ^Amaxln(cos(v3i< 1 (a:)))^ . 



The maxima (minima of the free energy density) are achieved at ot\ = ■ ■ ■ = = irn. 

In the absence of a quark mass, only the anomaly contribution in the imaginary part, 3F, of the free energy 
of a single domain appears under the logarithm of the cosine and defines the minima of the free energy density, 
ln(cos(w3i ;l (Q!))) = 0. Thus for massless quarks there is no continuous Ua(^) symmetry in the ground states, rather a 
discrete Z2 chiral symmetry. The anomaly plays a peculiar role here: selecting out those chiral angles which minimise 
the free energy so that the full Ua{^) group is no longer reflected in the vacuum degeneracy. It should be stressed 
here that this residual discrete degeneracy is sufficient to ensure a zero value for the quark condensate in the absence 
of a mass term or some other external chirality violating sources. 

Now switching on the quark mass, we see this discrete symmetry spontaneously broken, and one of the two vacua 
selected in the infinite volume limit according to the sign of the mass. In this case (for these conventions of boundary 
condition and mass term), the minimum at a = is selected. The quark condensate can be now extracted from the 
free energy via 

(i>(xU(x)) = - lim lim (vN)- 1 A e -»^("»). 

m^O N^oo dm 

Taking the thermodynamic limit N — > 00 and then m — > gives a nonzero condensate 

(i}(x)ip(x)) = < V>V >~ • (34) 

According to Eq. (31) the condensate is equal to 

(ip(x)tp(x)) = -(237.8 MeV) 3 (35) 

for the values of field strength B and domain radius R fixed earlier by consideration of the pure gluonic characteristics 
of the vacuum - string tension, topological succeptibility and gluon condensate. A nonzero condensate is generated 
without a continuous degeneracy of the ground states of the system. 

Dependence of the condensate on the domain radius R is illustrated in Fig.l. As expected the condensate diverges 
with R — > 00, since in this limit the number of low- lying strongly chiral modes is growing. We will return to the 
discussion of this feature again in the section on the meson spectrum. 



B. Multiflavour case: spontaneous breakdown of SUi,(Nf) x SUn(Nf) 

The question remains whether any continuous directions in the space of vacua are to be expected when the full 
flavour chiral symmetry is brought into play. For this we must generalise the analysis. We consider Nf massless quark 
flavours. Firstly, we observe that the fermion boundary condition in Eq.(5) explicitly breaks all chiral symmetries, 
flavour singlet and non-singlet (see also [7]). Thus the procedure we have used here of integrating over all a does 
not suffice to restore the full chiral symmetry of the massless QCD action. Rather, the boundary condition must be 
generalised to include flavour non-singlet angles, 

a -» a + f3 a T a , 

with T a the Nj — 1 generators of SU(Nf). Then integration over Nj angles a,(3 a £ [0,2n] must be performed for a 
fully chiral symmetric ensemble. The spectrum of the Dirac problem now proceeds quite analogously, except that the 
boundary condition mixes flavour components, thus an additional projection into flavour sectors is required in order 
to extract the eigenvalue equation analogous to Eq.(ll). 
For Nf = 2 the boundary condition can be chosen as 



(36) 
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1 2 3 

Domain radius R/R n 



FIG. 1: Absolute value of quark condensate as a function of domain radius (Ro = (760MeV) 1 ). 



where flavours now mix on the boundary. We must now solve the Dirac spectrum in the presence of this mixing. The 
spinors \ an d <P used in our solutions become now isospin doublets. We need to project the boundary condition onto 
separate equations for eigenvalues. 

Things proceed much as before in the spin sector with our previous decompositions and projectors. The boundary 
condition will devolve to the structure 



XW = e 



i(a+f3-3/2) 



<fi[\] 



where the 75 in the exponent is eliminated via the projection into chirality eigenspinors. The equation still mixes the 
flavour components and cannot be solved for an eigenvalue A. 
We introduce projectors 



P±(0) = 



l±0-5 



which have the property 



P±(P)e- i ^^ 2 = e ±i ^ 2 P ± (l3). 



Note here the appearance of the magnitude of the triplet [3 in the exponent. Projecting the spinors thereby we obtain 
separate equations for eigenvalues X[/3, ±] where we suppress all other quantum numbers previously dealt with, 

X [A[/3,±]] = e^ ± l^l/ 2 )^[A[/3,±]] 
It is convenient to define two new angles, 

a± = a± |/3|/2 

and in terms of these angles wc can calculate the quark determinant separately for each flavour projection as before. 
The free energy for each flavour is now a function of four constrained angles 

F± =F±{a,0 + ,0o,l3-) 

where we use this suggestive notation instead of /3i , /?2 , /?3 , and where ± denotes the isospin projection for the two 
quark flavours. For a given self-dual domain the total free energy will be a sum 

F T (aJ) =F{a+)+F{a-) 
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where F is the result from the one-flavour case. Under the assumption that the anomaly provides the entire chiral 
angular dependence of the determinant, the nonsinglet angles (3 a drop out due to the form 2q arctan(tan(a)) = 
2q(a + 2nir). This expresses the known result that the anomaly only depends on the flavour singlet directions or is 
Abclian. Thus for an ensemble of domains the free energy is identical to that for one massless flavour, namely it 
depends only on the Abelian angle a. Thus for Nf = 2, the Ua{^) direction remains fixed by energy minimisation 
while the SU(2)l x SU(2)r directions represent degeneracies in the space of ground states in the thermodynamic 
limit. 

For Nf = 3 the boundary condition will now involve the flavour mixing matrix 

e i/3 A a /2 

in terms of the Gcll-Mann matrices. It is now harder to explicitly diagonalise and project into flavour sectors. 
Nonetheless the form of the result is clear. The Cartan subalgebra consists of the diagonal generators A 3 and A 8 . 
Thus the result of a diagonalisation of the argument of the exponential will have the form, 

p a X a /2 -» b A \ A /2, A = 3,8 

where b l = b l (/3 a ), functions of the original nonsinglet angles analogous to |/3| for Nf = 2. Thus diagonalisation will 
amount to, 



e^^/ 2 ->diag( e b3+ ^ fc8 , e - b3+ ^ b8 , e -* fcS 



So projection of the combined U(l) x SU(3) chiral boundary condition will lead to three sets of equations with the 
eigenvalues depending on the combinations of angles 

B l = a + & 3 + -L& 8 , 
v3 

B 2 = a -b 3 + -^b 8 , 
v3 

B ° = °-7S" S 



where evidently 



^ Bi = 3a 



reflecting ultimately the tracelessness of the generators. The free energy for a given domain will be the sum of three 
terms 

3 

F T = ]TF(i? 4 ) 

i=i 

with F being the same expression evaluated for the one-flavour case. Once again, the assumption of the anomaly 
means cancellation of the (3 a dependent functions b A leaving only the a dependence in the free energy. Once again 
energy minimisation constraints the 1/4(1) direction, leaving the SU(3)l x SU(3)r directions unconstrained. Thus 
for Nf — 3 one expects eight, not nine, continuous directions in the space of vacua. 

The argument for arbitrary Nf is now evident, the key ingredients being the tracelessness of SU(Nf) generators 
and the specific form of the anomaly for the quark determinant. 



VI. ESTIMATION OF MESON MASSES. 



The consequences of the above realisation of chiral symmetries should be seen in the meson spectrum providing the 
splitting between pseudoscalar and vector mesons and between the pseudoscalar octet and 77'. 

The source for the difference between the masses of the octet states and -q' can be recognised in the drastic difference 
between correlators of the flavour octet Jp(x) and singlet Jp(x) pseudoscalar quark currents as they appear in the 
domain model, 



(J£(x)4(y)) = ((J£(x)4(y))) 

(Mx)My)) = ((Mx)Mv)))-((Mx)))((Mv)))- 



(37) 
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Here double brackets denote integration over quantum fluctuation fields and the overline means integration over all 
configurations in the domain ensemble. The second term in the RHS of the flavour singlet correlator contains two 
quark loops and is subleading in l/N c compared with the first one-loop contribution. The second term is entirely 
determined by the correlation function of the background gluon field B in the ensemble (8) and is proportional to the 
quark condensate squared. The analogous two-loop term in the flavour nonsinglet correlator is equal to zero due to 
the trace over flavour indices. It should be added that the pseudoscalar condensate, (tp^ip), naturally vanishes since 
parity is not broken in the ensemble of domains. Thus massless modes can be expected in the nonsinglet channel, but 
not in the flavour singlet due to the additional term in the correlator. This general structure of correlators is exactly 
the same as in the instanton liquid model [15] and manifests the mechanism for eta-prime mass generation proposed 
by Witten in [11] and appreciated in chiral perturbation theory by [16]. The correspondence between domain model 
correlators to original QCD diagrams is illustrated in Fig. 2. The grey background denotes averaging over the domain 
ensemble and is intended to represent the nonperturbative intemediate range part of the gluon exchange in the original 
QCD diagrams (a') and (b'). Diagrams (a) and (b) correspond to the first and second terms in the RHS of Eq. (37). 
Averaging in diagram (b) relates to both quark loops, which is the domain model representation of the exchange 
between quark loops in (b') by infinitely many gluons in the original QCD representation. It should be noted that 
diagrams (a) and (b) in Fig. 2 are the lowest order contributions in the fluctuation gluon fields Q which are treated 
as perturbations of the background. Higher orders include exchange by gluonic fluctuations in the presence of the 
domain mean field. 




FIG. 2: Diagrammatic representation of two types of contributions to the flavour singlet correlator: (a) and (b) as in the 
domain model; (a') and (b') as in full QCD. The grey background in (a,b) represents the nonperturbative part of the gluon 
exchange indicated in diagrams (a',b'). 

Direct calculation of meson masses within the domain model is not available yet. However we can estimate the 
effect of contributions subleading in 1 /N c in the flavour singlet pseudoscalar correlator within a calculational scheme 
which is quite close to the domain model in this respect. This is the main purpose of this section. Simultaneously we 
will schematically expose the form of the effective action for collective colourless meson-like modes as is expected to 
emerge in the domain model and estimate the value of the typical dimensionlcss parameter BR 2 but now from the 
meson spectrum. The model which will be used for this purpose is based on the bosonisation of a one gluon exchange 
interaction between quark currents and in which both quark and gluon propagators are exact solutions in the presence 
of a background (anti-)self-dual homogeneous gluon field as in the papers [17-19], but here with additonal correlations 
of the type (b) in Fig. 2. In [17-19], as well as the formulation of the model based on a homogeneous gluon field, 
applications to the calculation of the spectrum of light mesons, their orbital excitations, heavy-light mesons, heavy 
quarkonia, decay constants and form-factors are given. 

In the presence of a homogeneous (anti-)self-dual background field the quark condensate emerges due to chiral zero 
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modes, but, since there is a continuum of such modes in the homogeneous field, "an overkill" occurs [20, 21]: 

mfyijj) oc — B 2 , m — > 0, 

and the limit m — > cannot be defined properly. However, due to the same zero modes the momentum representation 
correlators of scalar, pseudoscalar, vector and axial vector currents have the following limits for m — > 

m 

flv/A^TBF 2 (p 2 /B). (38) 

Such a qualitatively different behaviour of correlators becomes manifest already at m 2 /B ~ 0(1) and leads to a strong 
splitting of the masses of the corresponding mesons ensuring light pions and heavy p-mesons in particular [18] (see 
also [15]). The fitted values of quark masses in the model correspond to the constituent masses. 

The definition of the chiral limit in that model thus cannot be given in terms of quark masses. However, as we 
discuss below it is nonetheless possible to define a regime when the calculated pion mass vanishes. The curious 
side of this approach to the description of chiral properties of the light mesons is that the zero-mode mechanism of 
condensate formation alone is exploited without any use made of a Nambu-Jona-Lasinio or, more generally, Dyson- 
Schwinger equation type mechanism [22]. The inclusion in this model of the additional correlations (b,b') of Fig. 2 
for the flavour singlet channel enables this approach to also work phenomenologically for r\ — rf masses with the chiral 
limit so-defined. 

In the domain model studied in the first part of this paper the quark condensate appears also due to specific 
chiral properties of (now nonzero) Dirac modes in gluon (domain-like) background fields and not due to a four- 
fcrmion interaction. The chiral properties of these low-lying modes resemble the properties of zero modes in the 
homogeneous-field model. The condensate diverges as R — > 00 as shown in Fig. 1 much like the condensate in the 
homogeneous field model diverges as 1/m. We expect that this dependence on R in the domain model should have 
a similar effect on meson correlators and masses as the 1/m singularity does in the homogeneous field case. The 
characteristic scale of higher modes is defined by the field strength B rather than the domain radius for BR 2 3> 1. 
Thus one would expect that the heuristic consideration below should be more consistent with the domain model 
picture if BR 2 > 1. 

The difference is that in the domain model flavour chiral symmetry is broken spontaneously, which is not the case 
in the homogeneous field model. So we should expect a massless pion for massless quarks: the current mass and 
the chiral limit are well-defined, as has been discussed in previous sections. However a detailed description of the 
emergence of Goldstone modes is a question of detailed study of correlators and the bound state problem in the 
domain model which is still beyond our efforts. 

To conclude these preliminary comments, the direct purpose of the following calculations is to demonstrate the 
(n, K) — r) — r}' splitting due to correlations Fig. 2b. An auxiliary purpose, justified by the similarity of chiral properties 
of zero modes in the homogeneous model and nonzero eigenmodes in the domain model, is a demonstration that the 
pseudoscalar singlet-octet mass splitting occurs simultaneously with a correct description of pseudoscalar- vector meson 
splitting. 




FIG. 3: Diagrammatic representation of the three-point vertex function Fq[q 2 q 3 - Diagram (a) corresponds to the first term 
in the RHS of Eq. (44), (b) and (c) to the second and third terms respectively. Solid lines are quark propagators S in the 
background field, vertices correspond to nonlocal quark-meson vertices Va i and the dashed lines represent meson fields $q 4 . 
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A. Effective action for colourless composite fields. 



We shall estimate meson masses within the model described by the following partition function 



l = 



g 2 Cgjf,(2) 



D 



r5e(-^|B), 



h Q 2 = ^W% 2 - 



-Mi 



(39) 

(40) 
(41) 



The effective action in Eq. (39) is expressed in terms of colourless composite meson fields &q(x) with the mass Mq 
defined by Eq. (40), where the condensed index Q denotes isotopic and space-time indices as well as all possible 
mesonic quantum numbers (iso-spin, spin-parity in the ground state, total momentum, radial quantum number), and 

fc-point nonlocal vertices Tg] Qk 



W fc [$] 



£ h Ql ...h Qk J d Xl ... J dx k ^ Ql {x 1 )...^ Qk {x k )T ( Q ) i Qk { Xl ,...,x k \B), 



Qi-Q 



G 



(i) 
Qi' 



r (3) 
C1S2S3 



T Qi Q 2 Q 3 S4 = G Qi Q 2 Q 3 qAxi,x 2 ,x s ,x 4 )-^E 2 {x 1 -X2) G%[ {x 1 )G%> 2 Q3 Q4 (x 2 ,x 3 ,x 4 ) 



G q1q 2 Q3( Xi > X2 > x ^ - o-2(xi - x 3 )G ( qI Q2 (x 1i x 2 )G ( qI{x 3 ) 



+ ^S 3 ( Xl ,X2, x 3 )GW { X1 )G%1 (x 2 )G^> (x 3 ) 



(42) 
(43) 

(44) 



(1)/ 



.(3) 



^ 2 (xi - a;3)G ( J i ) S2 (a;i,a;2)G^ Q4 (x 3 ,x 4 ) 



(2) 



+ E 3 ( Xl ,x 2 , x 3 )G%\ ( Xl )G^> (x 2 )G%> 3 Qa (x 3 , : 



--(2) 



- ^E 4 (x!,x 2 , x 3 ,x 4 )G ( qI {xi)G ( qI {x 2 )G [ q> 3 {x 3 )G [ q> 4 (x 4 ) 



2(1), 



,(1), 



(45) 



and analogous expressions for the higher vertices. Defining the meson-quark coupling constants Hq by Eq. (41) 
provides for the correct residue of the meson propagators at the poles and is known as a compositeness condition [23] . 
The vertices r( fc ) are expressed via quark loops Gq ^ with n quark-meson vertices 

G%l^ Qk ( Xl ,..., Xk ) = J da j TrV Ql (x 1 \B^)S(x 1 ,x 2 \B^)...V Qk (x k \B^)S(x k ,x 1 \B^) 

£ 

G Qi...Qi( Xl ' ■ ■ ■ > x G Q > l+1 ...Q k ( Xl + 1 > ■■■,Xk) = J da j 

s 

xTr{v Ql (x 1 \B^)S(x 1 ,X2\B^)...V Qk (x l \B^)S(xi,x 1 \B^)} 

x Tr [v Ql+1 (x l+1 \B^)S(x l+1 , x l+2 \BW )...V Qk (x k \B^)S(x k , x l+1 \B& ) } , (46) 
where bar denotes integration over all configurations of the background field with measure doj . The quark propagator 



S(x,y) = cxp ( --XpBpvyv J H(x-y), 



1 



p 2 /2vA 2 



1 - S 



m 2 /4vA 2 



1 + s 2 
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is the exact Dirac propagator in the presence of the (anti-)self-dual homogeneous field 

B^x) = - \^iB ilv x v , B^B^p = 4u 2 A 4 <5„ (9 , 

fl 

/a/3 = -^2 B ^' v = diag(l/6, 1/6, 1/3), 
where we introduced a scale A related to the field strength B as 

A* = — £?. 
2 

For arbitrary orbital momentum I and radial quantum number n the vertex function (more precisely - vertex 
differential operator) is factorised into a radial part 

l 

F nl (s) = J dtt l+n e st , s=V 2 /A 2 , 
o 

V//'= V -£/' V, tf = m f /(m f + m fl ), 

V M =a„ +iB„, V M = c^ -*-B„, (47) 

and angular part. In our particular calculation below we will not deal with the excited states. For the explicit form of 
the angular part of the vertex the reader is referred to [18], where technical details of the derivation of meson-quark 
vertices VQ k (xk\B^) and discussion of approximations and assumptions behind the derivation of the effective action 
for composite fields can be found. A simplified scalar field model allowing exact implementation of the method has 
been considered, and a variational procedure for approximate solution of realistic problems has been formulated in [24]. 
The relation of this method to the Bcthe-Salpeter equation and the nonrelativistic limit are analysed in [24, 25]. 

Thus the n-point vertex includes contributions of two types - the usual one-loop contributions (averaged over 
the background field) as diagram (a) in Fig. 3 and products of two or more such one-loop terms simultaneously 
averaged over the background field and multiplied by the corresponding correlators S of the background field, for 
example as in diagrams (b) and (c) of Fig. 3 for the three-point vertex. In the absence of a background field these 
additional terms would be just disconnected diagrams and would not appear in the effective action. For the purely 
constant field B they would correspond to infinite length correlations breaking the cluster decomposition property in 
the effective action (S n = 1). The idea of domains is implemented in the above expressions by means of finite length 
correlations of the background field which ensures cluster decomposition. 

The n-point correlation functions S„ of the background field are defined in Eq. (8). In particular the two-point 
correlator which will be used below can be written in explicit form 



N f 

S 2 (a;-j/) = — dz9(x - z)9(y - z) = 



2 +f( x -y) 



3tt V 4i? 2 



<Mp 2 ) 



3^ : 

— 3 arcsin(p) — 3pyl — p 2 — 2p(l — p 2 )\Jl — p 2 



(48) 



Geometrically, the correlation function — y) is equal to the volume of overlap between spherically symmetric 
four-dimensional regions determined by the characteristic function 9 and central points x, y, normalized to the volume 
of a single such region. 

Now all the elements of the effective action Eq. (39) are fixed: nonlocal meson-quark vertices Vq 1 (xi\B^) and 
quark propagators S(xk, Xi\B^) and background field correlators arc given in explicit analytical form. 

It should be noted here that meson-quark vertices are determined by the gluon propagator in the presence of the 
background field. In momentum representation both quark propagator, meson-quark vertices and correlators of the 
background field are entire functions in the complex momentum plane, which means that quarks and gluons are 
absent in the model as asymptotic states, which is treated as colour confinement. In the momentum representation 
the vertices have no imaginary parts and the probability of a meson to decay into a quark-antiquark pair is zero. 
As is shown in [17, 18] within the model formulated above and clarified in a simplified exactly solvable model in [24] 
entire propagators and vertices lead to a Regge spectrum of excited bound states represented by composite fields $g. 

The only free parameters of the model are the same as in the domain model discussed in the first part of this paper: 
the background field strength B, mean domain radius R, strong coupling constant g and quark masses mj. The factor 
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Cq in front of the gauge coupling constant g in Eqs.(39) and (40) is known explicitly 

Cjnl = Cj xl\{lln)V Cs ' p = \' Cs ' p = iV 

The only place where the gauge coupling constant g enters this scheme is the equation for meson masses Eq. (40); 
the rest of the effective action contains the meson-quark coupling constants calculated by means of Eq. (41). 

In the next subsection we present results for the masses of pseudoscalar and vector nonets coming straightforwardly 
from this formalism with a special emphasis on the r\ and i]' . 



B. The masses and decay constants of r) and r\ . 

With this representation, the calculation of masses of light mesons and weak decay constants, as performed in [18], 
is modified only with respect to the r\ and 7/ masses. Equations for these masses contain terms with additional 
correlation S 2 , and this splits them from the remaining light pseudoscalar mesons. 

As follows from Eqs. (40), the second term of (42) and (46), the simplest quark loops relevant to the calculation of 
masses are the constants 

G$ = -Tr\ a h 5 F 00 (x\B^)S(x,x\B^), (49) 

corresponding to the quark loops in Fig. 2b. Here a is the flavour index (a = 0, ... ,8). In the iso-singlet case (a = 0) 
this quantity is nonzero. If the s-quark mass is not degenerate with the masses of u— and d-quarks then Ggp is also 
nonzero. For a = 1, . . . , 7 the constants are identicaly equal to zero. 

Calculation of G^p is quite simple. The result of the action of the vertex operator F 00 (x\B) onto the quark 
propagator can be found as follows (in our case £/ = 1/2 in Eq. (47)): 

1 

F 00 {x\B)S(x,x\B) = y*dt expjJL (y x - Wy)^ S(y,x\B)\ x=y 




1 





Furthermore, 

1 



™> J (i^) ^ <*) 



Here the sign "=f" relates to self- and anti-self-dual configurations of the vacuum field. 

Substitution of Eqs. (51) and (50) in Eq. (49) and integration over the loop momentum p and the parameter t leads 
to the following integral representation 

t v 

where summation over / and v corresponds to traces in the flavour and colour spaces. As indicated in Eq. (47), v is 
a diagonal matrix relating to the direction of the vacuum field in colour space. 

Eq. (52) results in the following momentum representation of the two-point correlation function, Eq. (42) 

^% P (p 2 ) = n ab ( P 2 ) = u ab ( P 2 ) - su ab ( P 2 ), 
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where 



sn ab ( P 2 ) = -—{KRfT ab F{ P 2 ), T ab = x a f f x rr R f R f< 



ff 



F(p 2 ) = J dty/l - t 2 J dss cos (Wfl 2 t 2 s 



3tt 



3 arcsin \fs — (5 — 2s) \J s(l — s) 



T m — -(2R U + R s ) 2 , T 88 — — ^s) 2 ' — + R s ){Ru — R s ), 

n 00 = i(2n u + n s ), n 88 = ^(n„ + 2n s ), n 08 = ^(n„ - n s ). 

The function n^(p 2 ) has been calculated in Ref. [18], 

/ 

»^'- - W *i *j[ ds 'i te (^£T (its 



— X 

1 -i ■ A - 



P 2 Fl(tl,t2,*l,«2) . m / 



+ 



^2(S1,S 2 ) 



A 2 ipj(tl,t2,Sl,S 2 ) A 2 (1 - S 2 )(l - S 2 ,)<fil(tl,t2,Sl,S 2 ) 



+ 



2v{l-4vH 1 t 2 )F 3 {s 1 ,s 2 ) 



(fi 2 (h,t 2 ,s 1 ,s 2 ) 



cxp< - 



2wA 2 



<fi{h,t 2 , si,s 2 ) 



<Pi(ti,t 2 ,s 1 ,s 2 ) 1 , u , 

¥Wl,*2,Sl,S 2 ) 2 

¥? 2 = 2u(t 1 +t 2 )(l + sis 2 ) + (l + 4u 2 t 1 t 2 )(s 1 +s 2 ), 

= (1 + s 1 s 2 )[(s 1 +v(t 1 +t 2 ))(s 2 +v(t 1 +t 2 )) + v 2 (t 1 -t 2 ) 2 s 1 s 2 }, 
F 2 = (l + s lS2 ) 2 , F 3 = 2(l + sis 2 ). 

The masses of 77 and ?/ can be calculated by means of Eqs. (40), which take the form 



(53) 



1 + |-IM-M 2 ) = 0, 1 + |-IV(-M 2 ) = 0, 



^ 2 



n o + n 88 - v / (noo-n 8 8) 2 + 4ii 2 8 



n 00 + n 88 + yj (n 00 - n 88 ) 2 + 4n 2 8 



where II = II— 511 includes both diagrams in Fig. 2. These equations can be written more transparently if we introduce 
the mixing angle 



7] = r]s cos 9 + T/o sin 9, 
f\ = r]o cos 9 — 7]% sin 9, 

where the angle 9 is a function of momentum, 

tan2% 2 ) = 2n 08 (p 2 )/ (n 88 (p 2 ) - IIoo(p 2 

The polarization functions of r\ and rj and their derivatives with respect to p 2 take the form 

II,, = n 88 cos 2 9 + n 00 sin 2 9 + n 08 sin 26, 
U v , = n 00 cos 2 9 + fi 88 sin 2 9 - fl 08 sin 29, 
n' n = n 88 cos 2 9 + n' 00 sin 2 9 + n' 08 sin 29, 
n' v , = fl' 00 cos 2 9 + fl' 88 sin 2 9 - Yl' 08 sin 20, 

where for computing derivatives Eq. (55) has been used. 



(54) 



(55) 



(56) 
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With the values of the parameters 

m u = m d = 177.85 MeV, m s = 400.98 MeV, Vb = 469.52 MeV, g = 8.94, 

fitted from the masses of n, p, K and K* mesons such that 

= 140MeV, M K = 496MeV, U = 129.9MeV, f K = 150.8MeV, 
M p = M w = 770MeV, M K , = 890MeV, M = 1035McV, 



(57) 



(58) 



we find that 



for 



M v = 640 MeV, M v , = 950 MeV, 
h v = 4.72, h v > = 2.55 



VBR= 1.56. 



(59) 



The parameter R was fitted to rf mass. The domain model value for this dimensionless parameter quoted in the first 
part of the paper \J~BR = 1.24 is close to Eq.(59). Thus we see that all features of the spectrum of light vector and 
pseudoscalar mesons usually associated with chiral symmetries are correctly reproduced by the model quantitatively. 
The origin of the splitting between the pseudoscalar and vector mesons was discussed in detail in [18]. The new feature 
is the shift in the masses of rf and r\ with respect to the other pseudoscalar mesons and their mutual splitting, which 
occurs here due to the additional contribution of Eq. (53) to their polarization functions. The momentum dependent 
mixing angle takes different values at the scale of r\ and rf masses 



6 V = 0(-M%) = -19.8°, V , = 6{-M$) = 46.1°. 



(60) 



It is notable that 9 n coincides with the value of the mixing angle in the naive quark model, while V > is completely 
different from d n both in sign and magnitude. 

We quote also the result for the weak decay constants of 77 and rf : 

f° = 4.13 MeV, /8 = 165.7 MeV, = 154.5 MeV, 
f°, = 288.6 MeV, f%, = 23.67 MeV, f e q , = 183.3 MeV. 

These constants are defined by the matrix elements 



and /? corresponds to the mixings: 



(0\JZ li (0)\4>(p))=ip l *f$/V2, 

- Ql^t a q, t° = A°/2, t s = A 8 /2, 



Here 



/* = /S C os0 + /°sin0, 
$=$coB0-f$sm6. 

/» - hr, [fi cos + fS sin 0} , /«, - h v , If 3 cos - / 8 a sin 0} 
fo = g(2/ti + fs), f$ = g(/ti + 2/ s ), /§ = —^-(fu — fs)- 



The function //(— M 2 ) was calculated in [18], 



f m f sr 



4tt 2 



dtds\ds2{l + S1S2) 



2vt 



[si + s 2 + 2vt(l + sis 2 )] 3 
si s 2 



(l- Sl )(l- 32 ) 

(1 + Sl )(l+S 2 ) 



m 2 f /iv\ 2 



1 



+ 



1 



( M 2 T / 

CXP ( t,Sl,S2) 



* = 



sis 2 + vt(si + s 2 )/2 
si + s 2 + 2vt(l + S1S2) ' 
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As we have already mentioned, the quark masses here should be considered as constituent quark masses. The massless 
limit is ill-defined due to the contribution of zero modes to the quark propagator. However, as is discussed in [18], 
a peculiar feature is that all the necessary shifts and splittings in the meson spectrum occur explicitly due to the 
strong (and, for pseudoscalar, scalar,vector and axial vector states, very different) dependence of meson masses on the 
quark masses, which is driven entirely by the zero modes. In order to visualize this picture in a quantitative manner, 
including now also the 77 and 77', let us consider a regime in the model which can be called "a chiral limit". 

C. Massless pseudoscalar octet. 

This limit can be defined in terms of composite meson fields as the condition that the tt and K mesons become 
massless particles. Within the above-formulated meson theory this requirement is fulfilled if constituent quark masses 
satisfy the relation 

m s = rrid = m u = m* 

and the mass to* is defined from the equation for pion mass with m, = substituted, that is 

1 + ^ n ^ (0) = - (61) 

Since in this limit flavour SU(3) becomes exact symmetry the kaons are also massless. Furthermore, as follows from 
Eqs. (53) 

IIos = T 08 = T 88 ee 0, 

which means that the 77 meson is degenerate with the 7r and K. The mixing angle between 770 and 77s is equal to zero 
and the 77' meson corresponds to a pure flavour singlet state. As seen from Eqs. (53), 

<ffloo oc T 00 F(p 2 ) 

and the mass of the singlet state is split from the flavour octet states. Vector mesons are massive in this regime. 
Numerically one finds 

to* = 168.83 MeV, 

M v = M„ = M K = 0, M v , = 897 MeV, 

h n = h n = Iik = 4.28, h v > — .63, 

4 = fx = fn = 127.76 MeV, /° = 281.87 MeV, 

f v = # = 0. 

and 

M p = M K > = M u = 762MeV, M = 950MeV. 

This result displays a typical chiral limit picture with both flavour SUr(3) x 5/7l(3) and Ua{^) symmetries correctly 
implemented de facto: while the masses of octet states are drastically reduced from physical values to zero their decay 
constants are subject only to minor change. Simultaneously we see that the 77' does not look like a "Goldstone" 
particle: both its mass and decay constant /°, are practically unchanged. The vector states are subject to minor 
change also. For completeness we mention that the scalar and axial vector particles as ground state mesons are absent 
in the spectrum in this model, but they appear in the hypcrfinc structure of orbital excitations of vector mesons with 
quantitatively correct masses [18], which stays unchanged also in the chiral limit described here. 

In this picture m* looks like the condensate part of the constituent quark mass. This impression can be enhanced 
by estimating the current quark masses and their ratio, 

li u =m u -m* = 9.02 MeV, /_t s = m s - m* = 232.15 MeV, 
fjk s /fjk u = 25.73. 

The ratio, which is the only meaningful quantity, practically coincides with the generally accepted value fi s /^u — 25. 
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VII. CONCLUSIONS AND DISCUSSION 

We have shown that the correlation between the chirality of low- lying Dirac modes and the duality of the domain-like 
background field indeed drives the spontaneous breakdown of flavour chiral symmetry, as was hinted at in [4] , and as 
indicated in lattice calculations [5] . We have also shown that the mechanism generating the 77 — 77' mass splitting is the 
same as that causing area law confinement in the domain model, namely the finite range correlations induced by the 
domain mean field. In more detail, we have extracted the parity-odd part of the logarithm of the quark determinant 
and seen that the axial anomaly is recovered. The chirality properties of the Dirac modes in domains generate the 
anomaly in this context. We see however that in this formulation with domains with bag-like boundary conditions 
both the symmetric zeta and the asymmetric eta functions are necessary to obtain this result. We have then explored 
the consequences of the anomaly for the realisation of chiral symmetry in the domain model. The contribution of 
the anomaly to the free energy of an ensemble of domains leads to a spontaneous breaking of Z(2) symmetry rather 
than a £/a(1) for one massless fermion flavour. Quark condensation occurs because of the spectral asymmetry coming 
from the bag-like boundary conditions. In the multiflavour case, the Abelian nature of the axial anomaly guarantees 
that the discrete Z(2) symmetry remains spontaneously broken in addition to the correct continuous non-singlet 
SU(Nf) L x SU(Nf) R symmetries. One would thus naively expect only N 2 — 1 Goldstone bosons based on the 
number of continuous degenerate directions in the space of ground states of an ensemble of domains. 

On this basis we expect massless pions, kaons and eta mesons in the chiral limit. The domain model manifests 
correlations Fig. 2b in the singlet channel which generates the splitting in 77 — rf masses via the Witten mechanism. 
The contribution, Fig. 2b, is entirely driven by the correlators of the background field Eq. (48). Simultaneously the 
background field correlators are also entirely responsible for the area law [3] . Thus the origin of both mechanisms is 
identical. In the N c — > oo limit the rj' is massless since the contribution of Fig. 2b is 1/N C suppressed with respect to 
Fig.2a. 

We have seen that the calculational scheme does work well phenomenologically, but the dimensionless parameter 
turns out to be BR 2 w 0(1) such that for a self-consistent consideration of meson physics one needs to use the 
propagators of the domain background field: neglecting boundary conditions inside quark loops is not consistent with 
BR 2 rts 0(1). There is no clear separation of two scales as has already been observed in the consideration of static 
parameters of the vacuum - quark and gluon condensates, the string constant and the topological susceptibility within 
the domain model. 

Nevertheless, the calculations presented in the final section resemble two features of the domain model we have 
addressed in the first part of the article. The first feature is quite obvious: in both cases the additional contributions 
to the correlators (polarization functions) for 77' and, if m u ^ m s also for 77, are crucial for their splitting from the 
other pseudoscalar mesons, thus resolving this aspect of the 17a (1) problem. Certainly, the successful (but not self- 
consistent) quantitative description provided by the purely homogeneous background field need not ad hoc be equally 
successful in the domain model, and verification of this is one of our first priorities. The second feature is not quite so 
obvious. The splitting between pseudoscalar and vector meson masses in the case of the purely homogeneous field is 
determined by the singular behaviour of the quark condensate for m — > (as mentioned, the condensate diverges in 
the massless limit due to a continuum of zero modes [20, 21]). This singular behaviour is not present in the domain 
model - zero modes do not exist at all and the limit m — > is regular. However a quark condensate is generated 
in the domain model by the asymmetry in the spectrum of the Dirac operator, and this condensate diverges for 
R —* 00 as follows from Eqs. (31) and (34) and as shown in Fig. 1, as discussed above. This divergence is expected 
to play the same role for correlators of the domain model as m — > in the homogeneous field, thus generating strong 
pseudoscalar- vector splitting. The value of quark condensate given in Eq. (35) corresponds to R/Rq = 1 in Fig. 1. 
But, unlike the model based on a purely homogeneous field, the domain model manifests spontaneous breaking of 
flavour chiral symmetry and hence has the potential to reproduce a genuine picture of the chiral limit. This is only 
mimicked in the homogeneous field considerations above, though in a surprisingly successful manner. 

The zeta function calculation nonetheless remains incomplete, where the contributions that should eliminate the 
chiral angle dependent parity even part of the quark determinant's logarithm need to be calculated. The observation 
of lack of separation of scales, mentioned above, also means that a more careful study of collective modes within the 
domain model per se is necessary to put the phenomenological results above on a sound footing. In the same context, 
explicit correlation functions in the domain model would enable a study of the (anomalous) Ward identities and issues 
[26] related to the realisation of Goldstone's theorem (or otherwise) in the Ua(1) sector. 

There are two issues naturally related to this work and requiring at least some preliminary comment within the 
constraints of the calculations realised thus far. 

At first superficial glance, the mechanism of quark condensate generation through the asymmetry in the spectrum 
of the Dirac operator looks to be different to the approach of Banks and Casher [27] where the spectrum is symmetric 
and the condensate arises due to a finite density of chiral zero modes. However, these seemingly different formulations 
in fact have more in common. In particular, in the thermodynamic limit the system is characterised by a finite 
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density of low-lying (but non-zero) modes which are strongly chiral. Even if for a given fixed domain the spectrum 
is asymmetric, as indicated by the spectral asymmetry function 77, on average all positive and negative non-zero 
eigenvalues appear in the ensemble in a symmetric way, and neither left nor right chirality modes prevail in the 
ensemble. This can be seen in the symmetric distribution of the local chirality parameter given in our previous work 
[4]. The chiral condensate arises due to the existence of a finite density of strongly chiral eigenmodes both in the 
Banks-Casher formulation and in the model under consideration. The difference is that, in the domain model the 
corresponding eigenvalues are nonzero in the thermodynamic limit and none of the fermionic modes are purely chiral, 
rather they can be characterised by their average chirality with a definite sign correlated with the duality of the mean 
gluonic field in the domain: the lower the Dirac operator eigenmode the closer this average chirality to ±1. Thus, we 
would take the liberty to say that the domain model gives a "smeared" realisation of the Banks-Casher scenario. 

Another interesting question inspired by the chiral boundary condition used in the model is the manifestation or 
otherwise of features seen in chiral bag models of the nucleon. How much is in common here beyond the similar 
boundary condition? There is no simple answer to this question yet. We recall that the domain model is formulated 
in Euclidean four-dimensional space, colourless hadrons (if any) are not associated with domains themselves but are 
anticipated to arise as collective excitations of quantum fluctuations of quark (and possibly gluon) fields in the domain 
ensemble. The description of colourless hadrons requires analytical continuation to physical Minkowski space. Quark 
and gluon fluctuations in this approach are localised both in space and time: no asymptotic particle-like states can be 
associated with them. The model requires substantial use of methods of nonlocal quantum field theory. Various issues 
related to quantisation, unitarity, causality, Froissart-type bounds at high energies and interpetation of nonlocal fields 
as fluctuations localised in space and time can be found in [23, 24, 28-30] . The nonlocality appears here primarily due 
to the presence of strong background gluon fields which eliminate the pole in the momentum space quark propagator 
[17, 18, 31] rendering it an entire analytical function. The nature of coloured fluctuations as localised in space and time 
prohibits a straightforward resolving in the domain model of the manifestation of such interesting features as baryon 
number fractionalisation in chiral bags [32]. This phenomenon appears due to asymmetry (and a corresponding 77 
invariant) in the spectrum of the Dirac Hamiltonian, rather than the Euclidean Dirac operator, after a rearrangement 
of energy levels in the solitonic background chiral bag. There is no automatic one-to-one correspondence between 
phenomena in the domain model and those of the chiral bag model for the nucleon. In particular, the answer to 
the question of baryon number is intimately connected to the analytical properties of propagators (for example see 
[33]) and thus to the mode of dynamical confinement and the realisation of hadrons as propagating excitations in the 
domain ensemble, the subject of further work. 
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APPENDIX A: EXPANSION OF KUMMER FUNCTIONS 



1. Asymptotic form of Kummer function M '(k 2 : x 2 / '(4z) , k + 1, z) for k 3> 1, z, x — fixed. 



We use the representation [34] 



M(a/4 ^ + 1 ' z) = rR4^/ 




and change the integration variable = s giving 



klz 



4zofe+l-a/2z 



/ dss : 
Jo 



a/2z-k-l e -s*/4z Ik ( s y 



M(a/4z,k+l,z) 



r(a/4z) 



The asymptotic behaviour of 
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at a oc k 2 , k ^> 1 can be found by the saddle point method. Denoting a = a/2z, ( = 2z and f(s) = [1 — (k + 
l)/a] ln(s) — s 2 /2a(, we arrive at 



/>oo 

M fe = / dse af Wl k (s) 
Jo 



10 

with /(s) having a maximum at 

s = \/C(a - k - 1). 
Expanding the exponent about the saddle point and using the addition theorem 



z\ - z 2 e 



-i(f> 



zi - z 2 e 



to = \J 1 zf + z\ — 2z\z 2 cos((j)) 
with = and n one obtains the following representation 

-k\z~ °/4z2' c + 1— a / 2z 



M(a/z,k+ l,z) 



'2nz- 



[a - 2 Z (k + l)]«/^-(fe+ 1 )/ 2 e -a/4 Z+(fc +l)/2 e ,/4 



r(a/4z) 

CO 

x hn+k(so)In(z/^)[l + W 1 (x,z,n)/k + W 2 (x,z,n)/k 2 

71 — — CO 

+ £ J 2w+fc+ i(s )/ w (*/4) ^^ , 

n— — co 

z(z-2n + 3) 



(Al) 



W 1 = 



6x 2 



W 2 



72a; 4 



^(27z - 8xn + z 3 + \hz 2 + Azn 2 )^fl + x 2 + 8zn 2 x 2 - I2x 2 

+ I5z 2 - 9zx 2 + 2z 2 nx 2 + 8n 3 x 2 + z 3 x 2 + Az 2 x 2 - 56zn + 27 'z 
-A8n 2 x 2 + A6nx 2 + 26znx 2 + Azn 2 + z 3 ] 

bx 

s — \fa — 2z(k + 1), a = k 2 x 2 , x = tp, 

and calculation of the asymptotic form of M is reduced to the asymptotic decomposition of I 2n +k{so) an d ^2n+fc+i( s o) 
for k 3> n. 

Decomposition of the factors in front of the sum in Eq.(A2) at k 3> 1 gives 

J { \ z — a/4z2k+l~a/2z 



k\2 k 



{kx) k 



T(a/Az) 

^z/4+z/2x 2 
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The representation 



L(so) 
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/n2 v T(v + 1/2) 



dte- sot {i - t 2 y- 1 / 2 



is now suitable for determining an asymptotic decomposition of the I 2n +k(so) again by means of the saddle point 
approximation. With 
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and 



hi(x, n, z) 
h 2 {x, n, z) 



2{x 2 - 2)V1 + x 2 - 5 



24(1 + a; 2 ) 3 / 2 
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+ 17 - 144a; 2 + 276a; 4 + 52a; 6 + 384n - 576a; 2 n - 1152a; 4 n - 192a; 6 n - 576za; 2 - 576z] 
we obtain the following full expression for the asymptotics of the Kummer function up to 0(l/fc 2 ) 
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The series over n can now be resummcd in terms of the generating function of Bessel functions and its derivatives, 
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u = z/4, v 
Moreover, the presence of the factor 



1 + Va; 2 + 1 



, 1 = 0,1,2,3,4. 



,„ .„ 2 Va; 2 + 1 
exp < z/4 + z/2a; 2 - z— — - 2 



allows further simplification 
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The remaining formulae required can be obtained by differentiation. The final result is then 
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which reprodices the known asymptotic formulae [35]. This, as well as cancellation of powers of y in the denominators 
of Mi, gives a quite reliable criterion for correctness of the obtained expression. 

2. Asymptotic form of the Kummer function M(k 2 x 2 /4z + 1, k + 2, z) for k 3> 1, z, x fixed. 

Computations analogous to those of the previous section lead to the result 
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and 



M 2 (x,z) = —/ [385y 9 + (-840z + 2380)/ + (-5040z + 6048 + 840z 2 )y 7 

576(y +l) t> l 

+ (-12744z + 7636 - 480z 3 + 4480z 2 )y 6 

+ (9160z 2 + 144z 4 - 16896z - 1728z 3 + 3478)y 5 + (-3660 - 192z 4 + 6272z 2 - 10296z)y 4 

+ (-7440 - 32z 4 - 1192z 2 + 1728z 3 + 1008z)y 3 + (-1728z 2 - 5892 + 64z 4 + 480z 3 + 4680z)y 2 

+ (1728z + 16z 4 - 2439 + 408z 2 )y - 432 + 192z 2 ] 

with y = l/y/x 2 + 1. 



3. Asymptotic form of Kummer function M(—k 2 x 2 /4z, k + 1, —z) for k S> 1, z, x fixed. 

Use of previous calculations can be maximised by writing the function to be decomposed as 
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4. Asymptotic form of Kummer function M{— k 2 x 2 / (4z), k + 2, — z) for k 1, z fixed. 

The relation of M (-k 2 x 2 /(Az), k + 2, -z) to M(-k 2 x 2 / (Az),k+1, -z) is analogous to that of M(fc 2 x 2 /(4z) + 1, k + 
2, z) to M(k 2 x 2 1 (4z), k + l,z), thus the calculation is similar to that in Sec. A 3. The result is 

M(-a/4z,k + 2,-z) = e~ z M{k + 2 + a/Az, k + 2, z) 

k\2 k 1 



(kx) k V2^fc (1 + 

f Z yjx 2 + 1-1 



x exp • 



1 9 , 1 77 exp { k\J x 2 + 1 + k In ( f ] 1 

z 2 ) 1/4 I Vi + V^TT/i 

Mi(s,z) M 2 (x,z) 



where 



2 V^TT+i 



M (a;) = 



M (a:) + 



2y 



A; 2 



0(l/fc 3 ) 



1 + 
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M l {x,z) = - 



y(5y 5 + (12z + 22)y 4 + (48z + 12z 2 + 38)y 3 + (6 - 8z 2 - 



Ylz)y 1 + (-4z 2 - 48z - 39)y - 24) 



12(y + 1)3 



and 



M 2 (x,z) = 




j [385y 9 + (2380 + 840z)y 8 + (5040z + 6048 + 840z 2 )y 7 



+ 



+ 



+ 



+ 



(4480z 2 + 7636 + 12264z + 480z 3 )y 6 

(3478 + 1728z 3 + 8008z 2 + 144z 4 + 13152z)y 5 + (-2376z - 1792z 2 - 1152z 3 - 192z 4 - 3660)y 4 
(-19056z - 32z 4 - 11560z 2 - 2112z 3 - 7440)iy 3 + (-2880z 2 + 672z 3 - 13032z - 5892 + 64z 4 )y 2 
(384z 3 - 2439 + 864z + 2712z 2 + 16z 4 )y - 432 + 192z 2 + 2304z] 



with 



1 



y = 
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